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ABSTRACT 

Accurate photometric CoRoT space observations of a secondary seismological target, HD 174884, led to the discovery that this star 
is an astrophysically important double-lined eclipsing spectroscopic binary in an eccentric orbit (e ~ 0.3), unusual for its short 
3'?65705 orbital period. The high eccentricity, coupled with the orientation of the binary orbit in space, explains the very unusual 
observed light curve with strongly unequal primary and secondary eclipses having the depth ratio of 1-to-lOO in the CoRoT "seismo" 
passband. Without the high accuracy of the CoRoT photometry, the secondary eclipse, 1.5 mmag deep, would have gone unnoticed. A 
spectroscopic follow-up program provided 45 high dispersion spectra. The analysis of the CoRoT light curve was performed with an 
adapted version of PHOEBE that supports CoRoT passbands. The final solution was obtained by simultaneous fitting of the light and 
the radial velocity curves. Individual star spectra were derived by spectrum disentangling. The uncertainties of the fit were derived by 
bootstrap resampling and the solution uniqueness was tested by heuristic scanning. 

The results provide a consistent picture of the system composed of two late B stars. The Fourier analysis of the light curve fit residuals 
yields two components, with orbital frequency multiples and an amplitude of ~ 0. 1 mmag, which are tentatively interpreted as tidally 
induced pulsations. 

An extensive comparison with theoretical models is carried out by means of the Levenberg-Marquardt minimization technique and 
the discrepancy between models and the derived parameters is discussed. The best fitting models yield a young system age of 125 
million years which is consistent with the eccentric orbit and synchronous component rotation at periastron. 

Key words. Stars: binaries - Stars: variables: other - Stars: individual: HD 174884 



1. Introduction 



HD 174884 (Corot 7758) was selected as a seismology target 
for the first short run in the CoRoT "center" direction (SRcl) 
because - on the basis of its B8 V spectral classification - it 
was considered a good candidate to pinpoint the instability strip 
red e dge of the Slowly Pulsating B stars (e.g. IWaeUcens et al.l 
Il998l) . Its only available spectrum, obtained before launch with 



* Based on photometry collected by the CoRoT space mission and 
spectroscopy obtained with the CORALIE spectrograph attached to the 
1.2m Euler telescope at ESO, La Silla, Chile. The CoRoT space mission 
was developed and is operated by the French space agency CNES, with 
participation of ESA's RSSD and Science Programs, Austria, Belgium, 
Brazil, Germany and Spain 

** Postdoctoral fellow of the Fund for Scientific Research of Flanders 
(FWO) 



the Elodie spectrog raph and available from the Gaudi databas43 
(ISolano et al.l2005l) did not reveal the star's binary nature, which 
was instead immediately evident from CoRoT photometry. The 
system was continuously observed during 27 days. The light 
curve of HD 174884 is characterized by signatures of high ec- 
centricity: the secondary minimum markedly displaced from 
phase 0.5 and a bump, typical of close eccentric binaries, due 
to variable distortion of stellar surfaces along the orbit. The sys- 
tem is very interesting for several reasons: the high eccentricity 
is coupled with quite short a period (Porh = 3'?66), putting the 
system in an extreme location of the period - eccentricity dia- 
gram of B stars. This indicates young age or an inefficient cir- 
cularization process. Secondly, the huge difference between the 



' GAUDI (sdc.laeff.inta.es/gaudi/) is the data archive of the ground- 
based asteroseismology programme of the CoRoT mission. The 
GAUDI system is maintained at LAEFF, which is part of the Space 
Science Division of INTA. 
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Fig. 1. Upper panel: the complete CoRoT light curve of HD 174884 as normalized flux vs HJD. Lower panel: a blow-up, spanning 
approximately three periods, showing the tiny secondary minimum (first occurrence at reduced HJD of around 4218) 



primary and secondary minima, the latter detected only thanks to 
CoRoT accuracy, and the overall shape and accuracy of the light 
curve is a challenge for light curve fitting methods and a test of 
their adequacy to solve light curves of this unprecedented qual- 
ity. Finally, only a thorough analysis of HD 174884 could yield 
information of the possible pulsational properties of the com- 
ponents. In the following sections we present the results of the 
analysis of CoRoT photometry and of ground-based high disper- 
sion spectroscopy. A critical analysis of parameter uncertainties 
and of the uniqueness of the light and radial velocity curve solu- 
tion is an essential part of the study. 

The well constrained physical parameters of the binary com- 
ponents are then compared with evolutionary models computed 
with the CLES code (Code Liegeois d' Evolution Stellaire, 
IScuflaire et alj|2008bl) with the aim of constraining the age and 
the structural parameters (chemical composition, overshooting). 
Finally the light curve fit residuals, with variations of a few 10 
mag, are analyzed. 

2. Observations of HD 174884 and their analysis 

2.1. CoRoT photometry and extant data 

HD 174884 is a relatively bright star and several photomet- 
ric measurements could be found in databases and catalogs. A 
spectral classification as B8 V appears in the Simbad database 
(operated at CDS, Strasbourg, France) and was confirmed by 



the Gaudi classification. Available magnitudes and colors are 
collected in Table [1] the Stromgren photomet ry is from the 
Gaudi database, the infrared ones from 2MASS dSkrutskie et al.l 
i2006l) . These magnitudes and colors were dereddened accord- 
ing to|Moon (1985), assumin g a ratio of total to selective ab - 
sorption AJE(b -y) ^ 4.28 dCrawford & Mandwewalall 19761) . 
For the 2MASS colors t he respective ratios were taken from 
Rie ke & Lebofskvl (Il985h . The colors, once corrected for the 
rather heavy absorption, correspond to the expected values for 
the (primary) spectral type, and in particular rule out the pres- 
ence of a near IR excess, which was suggested by the observed 
values. The star was also observed by Hipparcos, but its parallax 
is of little use, being affected by a large error, n - 1.95 + 0.73 
mas, acc ording to the new reduction of the Hipparcos data 
(Ivan Leeuwen 2007.) . 

CoRoT observed HD 174884 (Corot ID: 7758) for 27 dayfl 
The almost uninterrupted time series extends over about seven 
orbital cycles with a time sampling of 32^ The full light curve 
has 60016 points, after removal of the measurements perturbed 
by South Atlantic Anomaly crossing (proton events). Each point 
is the mean made of 32 1-s exposures and has, therefore, an as- 
sociated r.m.s. deviation, which is typically of the order of 10"^ 
times the measured flux. 



- The data can be retrieved from the CoRoT public archive at 
http://idoc-corot.ias.u-psud.fr/, 
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A slow decreasing trend, presumably due to some kind of in- 
strumental decay or gradual pointing drift, was corrected by fit- 
ting and subtracting a second order polynomial passing through 
the maxima. Furthermore, we deleted a few hundred spurious 
points that arise from the inadequate linear fit used in the first 
version of the CoRoT pipeline. The complete, normalized, light 
curve (59492 points) is shown in Fig. [T]together with a blow-up 
of the central part, displaying a tiny 1.5 mmag deep) sec- 
ondary minimum. 

The binary period was derived in two s teps: an initial es- 
timate with Period04 (iLenz & Bregerl |2005[) was refined wit h 
the Phase Dispersion Minimization method (Stellingwer fl 19781) . 
which is better suited for the analysis of non-sinusoidal varia- 
tions. The epoch of primary minimum was derived by a linear fit 
to the seven minimum epochs, which were obtained by parabolic 
fit to the curve minim a. The results were checked with PHOEBE 
dPrsa & Zwitted2005l) . Note that PHOEBE determines the epoch 
of zero phase with respect to the argument of periastron (rather 
than with the position of the primary minimum), and provides 
phases for important positions in orbit, like superior conjunc- 
tiorQ. We adopt the following ephemeris: 

r^ini = 2454204.68162(4) -i- 3'?65705(1) x E. (1) 

The light curve was folded with this period and found to be 
stable within the observed time frame, as the dispersion of the 
phased curve was comparable to that of the original time series. 
Therefore, we binned the oversampled curve into 650 normal 
points, to facilitate the light curve solution. The binning was per- 
formed with a variable step, A0 - 0.002 outside the eclipses and 
A(p - 0.0005 during the eclipses to keep good sampling. The 
mean number of points contributing to normals is ~ 130 outside 
the eclipse and ~ 30 inside, the folded light curve is shown in 
Fig.E] 

2.1 .1 . Preliminary liglit curve analysis 

The light curve fit procedure went through several stages, also 
because the spectroscopic data were obtained only some months 
after the CoRoT run, thanks to a follow-up program organized 
after the discovery of binarity. Our final fit is a simultaneous so- 
lution of both light and radial velocity curves, but the first steps, 
based on photometry alone, were quite instructive. They were a 
test of what can be derived from a single light curve which, on 
the one side, is of unprecedented accuracy but on the other has 
features suggesting weak constraining of the solution. 

All light curve so lutions were performed with PHOEBE 
jPrsa & ZwitteJ l2005l) . which was expanded to include flux 

^ In the case of an elliptical orbit, the minimum light epoch can be 
slightly shifted from the superior conjunction; in our case we observe a 
shift of ~ 4 minutes. 



Table 1. Colors and magnitudes of HD 174884 





Observed 


Dereddened 


V 


7.998+0.002 


7.001 ±0.006 


b-y 


0.173+0.001 


-0.055±0.002 


mi 


0.042 ±0.001 


0.117±0.001 


Cl 


0.637+0.005 


0.594±0.005 


P 


2.701 ±0.011 




J 


7.577+ 0.024 


7.302±0.025 


H 


7.562± 0.023 


7.391 ± 0.024 


K 


7.521± 0.021 


7.412±0.021 
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Fig. 2. Upper panel: the phased and binned CoRoT light curve 
of HD 174884. Dots: normal points as described in the text. 
The solid line is the light curve fit according to the model of 
Sect. [3] which can be distinguished from normal points only dur- 
ing primary eclipse. Lower panel: the corresponding residuals 
(observed - computed normalized flux) 



computation in the CoRoT passbands for both seismology and 
exo-planet field. Specific passband intensities and limb dark- 
ening coefficients were c omputed for CoRoT passbands from 
ICastelU & Kuruczl (l2004 's NEWODF models spanning from 
3500 K to 50,000 K across the entire H-R diagram. The spec- 
tral energy distribution functions (SEDs) were synthesized with 
SPECTRUM32.75 jGrav & Corballvlll994 at a 1 A dispersion, 
on a 3000 A to 10,000 A wavelength range. The limb darken- 
ing tables were computed by synthesizing SEDs in 20 values of 
H = cos 6 from the center of the disk to the limb, and fitting 
the linear cosine law and non-linear log and square root laws by 
least squares. 

For the preliminary solutions, we profited from the scripter 
version of PHOEBE and used Nelder & Mead's Simplex min- 
imization for quick descent towards the minima and stan- 
dard diff'erential corrections when close to the solution. The 
value of primary effective temperature was fixed to that pro- 
vided by Gaudi, which is derived by means of the code 
TempL ogCfl and Stromgren photo metry. TempLogG uses the 
grids of iMoon & Dworetskvl (|1985[) in the Teff-log g para meter 
space, with the improvements by [Napiwotzki et alj dl993 ?). The 
resulting value is Tg^ i = 13140 K; we did not use the internal 
error from TempLogG but assumed an uncertainty of 1500 K, 
expressing the a ccuracy of T^ff determi nation for B-stars by dif- 
ferent methods jDe Ridder et al.ll2004t) . For the secondary star 
we actually started from a reff" value (3500 K) very far from the 
final one (wrongly) assuming, on the basis of the great difference 
of in-eclipse depths, a small, cool secondary star 

We computed, as starting value of a grid, pairs of eccentric- 
ity (e) and longitude of periastron {a>) providing the observed 
distance between maxima. The parameters adjusted in the min- 
imizations were inclination /, e, lo, Tsff.2, and surface potentials 
Q.1,2', the primary passband luminosity Li was separately com- 
puted, rather than adjusted, as this allows a smoother conver- 
gence to the minimum (Prsa & Zwitteij|2005l) . For limb darken- 
ing we adopted a square root law that employs two coefficients, 

http://www.phys.appstate.edu/spectrum/spectrum.html 
http ://w w w. univie. ac . at/asap/templogg/main. php 
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xiB and yio per star and per passband. The gravity darkening 
and albedo coefficients were kept fixed at their theoretical val- 
ues^ = 1.0 and A = 1.0 for radiative envelopes (vonZeipe| 
|1924|) . The mass ratio was set to values between 0.1 and 0.5. 
When the adjustment of these parameters did not provide further 
improvement, we attempted to adjust synchronicity parameters 
(ratios of rotational to orbital velocity), Fi;i. 

In spite of the very shallow secondary eclipse and the low in- 
clination, the solution we obtained is fairly close to the final one, 
implying that the single light curve is already constraining the 
solution (this was confirmed later on by the results of heuristic 
scan, presented in Section l373T i. In particular, we could exclude 
solutions with a cool secondary, in favour of a system config- 
uration yielding a grazing eclipse of a relatively hot secondary 
star. The very different depth of eclipses is actually due to the 
combination of high eccentricity and orbit orientation in space. 

It was interesting to check, a posteriori, that this single light 
curve solution gave the following deviations for the various pa- 
rameters with respect to the combined solution of photometry 
and radial velocity data: 6% for e and a>, 1% for /, 15% for Tgff 2, 
6-15% for the component radii (ri,r2), 12% on Li. We did not 
adjust the mass ratio, and actually the true value is outside the 
range we explored, but we had indeed a tendency of smaller val- 
ues of the cost function for the higher q values. Furthermore, the 
value of Fi was within 1% of that derived later on from line pro- 
file fitting. On the other hand the value of F2 was almost 4 times 
the spectroscopic one. 

In conclusion the light curve solution was well constrained 
for the main parameters and for primary rotational distortion. 
This has to be put in relation to the high eccentricity signatures 
which act as pinpoints on the parameter values. 

2.2. Spectroscopic follow-up 

Given that HD 174884 turned out to be a close binary, a spec- 
troscopic campaign was organised as soon as possible after the 
CoRoT space photometry had become available to us. We used 
the CORALIE echelle spectrograph attached to the 1.2m Euler 
telescope at La Silla in Chile. The star was monitored from 7 
to 16 June 2008, adopting an exposure time of 30 minutes. We 
used the online reduction package available for the CORALIE 
spectrograph and based on the method by Baranne et al. (1996). 
The process involves the usual steps of de-biasing, flat-fielding, 
background subtraction and wavelength calibration by means of 
measurements of a ThAr calibration lamp. All spectra were sub- 
sequently normalized to the continuum by fitting a cubic spline 
and barycentric corrections were computed for the times of mid- 
exposure. This led to barycentric wavelength-calibrated normal- 
ized spectra with a S/N ratio between 50 and 120 in the blue 
spectral region (the red parts suffered more from relatively poor 
atmospheric conditions during some nights), nspection of the 
spectra revealed that HD 174884 is a double-lined spectroscopic 
binary, evident line doubling could be seen in Hei and, more 
clearly in Mg 11 4481 A line (see Fig. [3]), identifying both com- 
ponents as stars of spectral type B. The finding of similar tem- 
peratures from the preliminary light curve solution was therefore 
fuUy confirmed. 

The set of 45 spectra was independently analyzed with dif- 
ferent methods: profile fitting of the abovementioned Mgii line 
and Fourier spectra disentangling. The radial velocity curves 
from profile fitting were used in the simultaneous solution dis- 
cussed in Section 13.11 Fourier disentangling provided the indi- 
vidual stellar spectra and allowed, as well, a cross-check with 



Table 2. Radial velocities of HD 174884 



BJD -245000 




Vmd.l 


V,-ad.2 




(km s-') 


(kms ') 


4625.6910668 


-94.5 


±1.6 


85.1±2.6 


4625.7333949 


-97.1 


±1.7 


87.6±2.6 


4625.7767647 


-97.5 


±1.5 


99.5 ±2.5 


4625.8241393 


-104.4±1.5 


99.3±2.5 


4625.8549044 


-107.6 


±1.5 


99.0±2.5 


4626.6823677 


-78.7 


±1.6 


69.8±2.7 


4626.7347772 


-74.8 


±1.6 


63.9±2.7 


4626.7602527 


-76.1 


±1.6 


60.4±2.7 


4626.7880894 


-74.4 


±1.6 


58.2±2.8 


4626.8353945 


-70.3 


±1.6 


53.4±2.8 


4627.6872211 


24.9 


±2.0 


-83.9±2.9 


4627.7463070 


31.6 


±2.0 


-90.4±2.8 


4627.7684857 


35.1 


±2.0 


-98.2±2.7 


4627.8089965 


36.6 


±2.0 


-107.4±2.5 


4627.8536741 


43.6 


±1.9 


-113.1±2.5 


4628.6470598 


77.1 


±1.7 


-162.4±2.5 


4628.6771883 


66.7 


±1.8 


-152.9±2.5 


4628.7144813 


58.6 


±1.8 


-137.3±2.5 


4628.8275062 


23.3 


±2.1 


-91.6±2.8 


4628.8613847 


15.2 


±2.2 


-80.6±2.9 


4629.6750711 


-108.0+1.5 


105.0±2.5 


4629.7286030 


-106.9 


±1.5 


108.1±2.5 


4629.7752019 


-106.6 


.±1.5 


107.4±2.5 


4630.6845270 


-54.8 


±2.0 


30.7±3.5 


4630.7230468 


-48.7 


±2.0 


23.9±3.5 


4630.7712660 


-44.7 


±2.0 


18.1±3.5 


4630.8161634 


-42.1 


±2.0 


13.7±3.5 


4630.8495673 


-32.2 


±2.0 


5.1±3.5 


4631.7076510 


71.8 


±1.7 


-159.7±2.5 


4631.7415409 


77.3 


±1.6 


-175.2±2.5 


4631.7768777 


83.5 


±1.7 


-176.7±2.5 


4631.8041587 


85.6 


±1.7 


-182.6±2.5 


4631.8305253 


89.7 


±1.8 


-187.0±2.5 


4632.6861425 


-32.9 


±2.0 


3.2±3.5 


4o3/. /41Dyjo 


-51.7 


±2.0 


n 1 Q c 
/3.0±3.5 


4632.7934837 


-62.5 


±2.0 


45.0±3.0 


4632.8199197 


-72.8 


±1.8 


48.0±3.0 


4632.8472122 


-75.2 


±1.6 


53.9±3.0 


4633.7128864 


-98.8 


±1.6 


92.8±2.5 


4633.7401557 


-95.9 


±1.6 


94.8±2.5 


4633.7913378 


-94.1 


±1.6 


87.2±2.5 


4633.8186997 


-91.9 


±1.6 


81.4±2.5 



the results of radial velocity analysis, that turned out to be very 
useful, given the peculiarity of the system configuration. 

2.2.1 . Radial velocity curves from profile fitting 

We extracted the region 4450-4500 A and fitted the Mg n 448 1 A 
line profile with a standard non-linear least squares fit proce- 
dure (IDE's CURVEFIT) using a user-supplied non-linear func- 
tion with an arbitrary number of parameters. The supplied func- 
tion was a theoretical profile obtained by convolving for each 
star a Gaussian (corresponding to CORALIE spectra resolu- 
tion) and a (dominant) rotationally broadened profile. In prin- 
ciple there were 8 adjustable parameters (position, width, depth 
of each profile and the values of continuum on the two sides). 
However, the iterative procedure did not converge in the case 
of heavily blended profiles. Therefore, we fitted first the un- 
blended profiles, adjusting the full set of 8 parameters, and fixed 
then the rotational velocities, sin /, to the average values ob- 
tained in this way for the remaining fits. These average values 
arev'rjsin/= 109. 3 + 1.5 kms"' and Vi-_2 sin/ = 59.7+2.7 kms"'. 
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Fig. 3. An example of the fit of the line profiles of Mg n 448 1 A . The black solid line is the fit with a function obtained by summing 
the individual profiles of the components (a convolution of a gaussian and rotational broadening). For clarity a gray line is also 
shown, which corresponds to a smoothing with an averaging boxcar of width comparable to the FWHM of the interstellar Na ii 
D line (the narrowest lines appearing in the spectra). Left panel; unblended lines at phase ip - 0.71, right panel: blended lines at 
(f - 0.45. The lower boxes show the fit residuals. 
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Fig. 4. The radial velocity curves for the two components of 
HD 174884 (dots: primary , triangles: secondary component). 
The phases are computed according to the ephemeris in Eq. 
[U velocity unit is kms The solid line is the fit described in 
Sectionirn 



The fits were performed on the original profiles and as well 
on the smoothed ones, with different boxcars (see Fig.|3]for ex- 
amples). The results were very similar even with resolution de- 
graded down to 0.025 A/pix, i.e. the FWHM of the narrow Nal 
D interstellar lines (the narrowest lines seen in the spectra). The 
values in Table |2] correspond to the fit of the unsmoothed pro- 
files. The fit uncertainties were derived with a bootstrap experi- 
ment (see Section [XZt . providing a more meaningful statistical 
estimate than formal errors of the fit. 

Only 42 out of 45 spectra yielded a determination of radial 
velocities, for the remaining three, too close to conjunctions, the 
fitting procedure did not converge. The resulting values are col- 
lected in Table |2] and the radial velocity curves are shown in 



Fig Sl The radial velocity curves were fitted with PHOEBE. This 
preliminary solution was not simultaneous, but we used an itera- 
tive procedure between light and radial velocity curve solutions, 
feeding the output of the iteration on one set of data as input 
to the other The purpose was to get a set of input parameters 
for the final combined solution presented in Section [3] Actually 
the parameters obtained from the preliminary light curve fit of 
Section [2.1. II were found to predict fairly well the radial veloc- 
ity curves, apart of course from the mass ratio. The values best 
fitting the radial velocity curves alone were: e - 0.292 + 0.003, 
w = 51.85° +0.5, flsin/ = 18.27 + 0.08 Rq, y = -20.7 + 0.3 
kms"', q - 0.674 ± 0.005 (formal errors). The corresponding 
radial velocity amplitude for the primary is K\ - 106.4 + 0.6 
kms^^ 

2.2.2. Fourier disentangling of spectra 

Sou thworth & ClausenI (l2007h carefully compared various pos- 
sible present-day techniques to derive the orbital elements of 
spectroscopic binaries and concluded that the best results are ob- 
tained via properly applied disentangling. They showed that it is 
necessary to map the dependence of the sum of squares of resid- 
uals in the semiamplitudes K\ and K2 (or the mass ratio K1/K2) 
since this sum usually has a number of local minima in the pa- 
rameter space. 

To obtain reliable orbital elements of HD 174884 based on 
all 45 spectra we applied the disentangling in the Fourier do- 
main developed by iHadraval ([1991 [199^ iBol . using the latest 
publicly available version of his Fortran program KOREL, dated 
December 2, 2004. 

In gen eral, our experience wi th KOREL fully confirms the 
findings bv lSouth worth & Clausen (.2007 ) and also that the result 
is sensitive to the choice of starting values of the elements and 
initial values of the simplex steps, i.e. that the simplex solution 
can easily be trapped in a number of various local minima. 

The spectra at our disposal have a large resolution and it 
would be very time-consuming to carry out disentangling for 



6 



C. Maceroni et al.: HD 174884: a strongly eccentric, short-period early-type binary system discovered by CoRoT 





58700 




58600 




58500 


o 




"O 


58400 


'co 




<D 


58300 


i 


58200 






o 


581 00 


Z) 




O" 


58000 






o 


57900 


Z) 


57800 




57700 




57600 



3 blue segments 



100 102 104 106 1 08 110 112 114 
Ki (km s-1) 



65000 

64000 

I 63000 

^ 62000 
o 

$ 61000 
I 60000 
E 59000 
58000 
57000 



3 blue segments 



0.6 0.7 0.8 

K1/K2 (km s-1) 



1 0660 
1 0650 
1 0640 
1 0630 
1 0620 
1 061 
1 0600 
1 0590 
1 0580 
1 0570 
1 0560 
1 0550 
1 0540 



He I 5876 with telluric lines 



100 102 1 04 1 06 U 
Ki (km s-1) 



110 112 114 




0.6 0.7 0.8 

K1/K2 (km s-1) 



Fig. 5. A search for optimal values of Ki and the mass ratio Ki jKi- 



the whole interval of wavelengths they cover. At the same time, 
the S IN of the spectra is not particularly high. To cope at least 
partly with this problem and to obtain good disentangled spec- 
tra containing spectral lines sensitive to the determination of at- 
mospheric paramaters and chemical abundances, we slightly de- 
graded the spectral resolution, averaging the spectra over 0.05 A 
and used only the following four spectral regions: 4050 - 4190, 
4450 - 4500, 48 10 - 4910, and 5860 - 5885 A. For each spectral 
segment, we estimated the S /N ratio of individual spectrograms 
and used the weights proportional to {S /N)^. The disentangling 
was then carried out separately for the three blue spectral re- 
gions combined into one solution and for the red region near 
the He I 5786 A line where we also considered the telluric lines 



Table 3. KOREL disentangling orbital solutions for the three 
combined blue spectral regions near some stronger spectral lines 
useful for determination of atmospheric parameters, and for the 
red spectral region near the He 1 5876 A line. The period, eccen- 
tricity and longitude of periastron were kept fixed at the values 
of 3'?65705, 0.29394, and 51?31, respectively. The epoch of the 
periastron passage Tperiast. is tabulated in HJD-2454628. 



Element 


3 blue regions 


He I 5876 


^periast. 


0.69173 


0.69037 


K,(km s-1) 


106.983 


106.956 


A:2(kms-') 


159.39 


161.64 


KilKi 


0.6712 


0.6617 



in the solution. In all cases, we derived also the rel ative line 
strengths for individual spectra during the solutions dHadraval 
1997). In all solutions, we kept the orbital period P, eccentricity 
e and the longitude of periastron lo fixed at the values of 3'?65705, 
0.29394, and 51?31, respectively, which were derived from the 
final light curve fit (see next section). The epoch of the periastron 
passage Tperiast. was chosen inside the time interval covered by 
the spectra, near HJD-2454628. 7, but was allowed to converge. 

We first investigated the variation of the sum of squares of 
the residuals as a function of the radial velocity semi-amplitude 
of the primary only, keeping all other elements but the periastron 
passage fixed. This dependence was investigated separately for 
the combined blue, and the red spectral segment and is shown 
in a graphical form in the left panels of Fig. |5] This led to the 
conclusion that the true semi-amplitude of the primary must be 
close to 107 kms"'. Fixing this value of the semi-amplitude, we 
then investigated the sum of squares of residuals as a function 
of the mass ratio, again separately for the blue and red spectral 
segments. This also rather uniquely restricts the plausible values 
of the mass ratio to 0.66-0.67 (cf. the right panels of Fig.|5]). 

The final disentangling was then derived starting with the 
best values found and allowing for the convergence of K\, 
K1/K2, and Tperiast.- We did a number of trial solutions for each 
spectral segment, varying slightly the initial values of the el- 
ements and starting simplex steps to find out solutions with 
the smallest sums of residuals. The results are summarized in 
Table |3] Their comparison provides some idea about the prob- 
able errors of the elements. We remind that for the region 
5860 - 5885 A, we also included the weak telluric lines into the 
solution. The comparison of results is useful since the program 
KOREL we used does not provide any estimate of the errors of 
the elements. 
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Fig. 6. Disentangled line profiles. 



We included KOREL results in PHOEBE as cross-check, 
following a somewhat unusual procedure. For each of the four 
KOREL solutions we extracted radial velocities of both compo- 
nents for the 45 times of observations and derived their average 
values over the two solutions for each B JD, giving 3 times higher 



weights to the blue RVs based on the separate spectral segments. 
These "mean theoretical" RVs were used to verify the results 
from PHOEBE that are based on profile fitting RV curves; we 
find excellent agreement between the two methods of RV ex- 
traction. 

Keeping the elements from the final PHOEBE solution, dis- 
cussed in the next section, fixed, we disentangled the binary 
spectra. They are shown in Fig.|6] 



3. Final analysis 

3.1. Simultaneous light and radial velocity curve fit 

To obtain a set of parameters that most accurately describes 
physical and geometric properties of HD 174884, we submitted 
photometric and RV data in the time domain to a simultaneous 
fitting process in PHOEBE. Simultaneous fitting imposes con- 
sistency of the solution across all available data-sets. 

The initial values of Teffj and Teff.a were adopted from the 
previous preliminary fits: 13, 140K and 11, 000 K, respectively. 
The values of gravity darkening coefficients and albedoes were 
set to their theoretically expected values for radiative envelopes: 
gi,2 - 1-0; Ai 2 - 1.0. The center-of-mass radial velocity Vy was 
derived from the RV fit and held constant. 

We used the RV derived value of a sin/ = 18.27 + 
0.08 R0(cf. 92.2.1b to constrain the solution to those values of 
a and / that preserve a sin /. Line profile analysis and spectrum 
disentangling (cf. 92.2.2b yield consistent values of vi sin / = 
109.3 + 1.5 km s"' and V2sin/ = 59.7 + 2.3 km s"'; these were 
used to constrain synchronicity parameters Fij = cOrot/Morb, i-C- 
the ratio of rotation to mean orbital angular velocity. Square- 
root limb darkening coefficients are interpolated on-the-fly from 
PHOEBE lookup tables, based on the current values of Tetj, 
log g, [Ml H] and Vrot for each star. As before, passband luminos- 
ity was computed rather than adjusted, to enhance convergence. 
The constraints are imposed during each iteration, projecting out 
only those solutions that satisfy the embedded physical require- 
ments. 

Parameters that were adjusted are: eccentricity e and argu- 
ment of periastron u, inclination /, surface potentials Qi and Q2, 
mass ratio q, and secondary surface temperature Teft-.a. The mass 
ratio typically has little detectable signature in well detached EB 
light curves since it affects ellipsoidal variations; in the case of 
HD 174884, however, the CoRoT light curve's sensitivity on q is 
significant, due to the eccentric orbit, the consequent star shape 
distortion and high data accuracy. These factors allowed us to 
adjust the value of q based on simultaneous LC and RV data. 
There is an excellent agreement between the q value from the fi- 
nal fit and that derived from KOREL. The difference is less than 
1-sigma interval for the three blue segments and within 2 sigma 
for the red one (which has a lower S/N ratio). 

The final parameter set is listed in Table |4] 

Before adopting our final solution and proceding to compari- 
son with the stellar models we felt there were still two questions 
to address: solution uniqueness and a realistic estimate of the 
uncertainties. 

It is well known that light curve solutions can be affected 
by uniqueness problems, and this could be a serious issue in the 
case of our light-curve that features a grazing eclipse. It is also 
well known that the formal errors derived from least squares co- 
variance matrix are inaccurate because the parameter space axes 
are not orthogonal and marginalization around the solution is 
not fully justified. More realistic errors were derived by boot- 
strap and by full-fledged heuristic scanning, performed on all 
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Table 4. Physical and geometric parameters of HD 174884. 



Parameter: 


Primary 


System 


Secondary 






18.88 ±0.14 




in 




75.35 ± 0.03 








0.674 ± 0.009 




e 




0.2939 ± 0.0005 




(X) 




51.31" ±0.06° 




Vy (km s"') 




-20.6 ± 0.3 






6.055 ± 0.008 




7.866 ± 0.020 


Teff (K) 


13140 




12044 ± 50 


^(Mo) 


4.04 ±0.11 




2.72 ±0.11 




3.77 ± 0.03 




2.04 ± 0.02 




-1.71 ±0.04 




0.01 ± 0.05 




2.0* ± 0.2* 




2.3* ±0.2* 



*computed from the observed v sin i and orbital period. 



adjusted parameters. The latter procedure provides information 
on the morphology of the cost function hyper-surface in the ad- 
justed parameter space. The final estimates of most parameter 
uncertainties appearing in Table |4] is obtained by the heuristic 
scan presented in Section 13.31 As, however, this is a procedure 
rather expensive in terms of computing time and resources , we 
found it instructive to compare the results of both methods, and 
quantify the differences. Furthermore, we adopted the uncertain- 
ties from bootstrap for the parameters not adjusted in the scan. 

3.2. Parameter uncertainties by bootstrap 

Bootstrap resampling is a very useful technique to estimate pa- 
rameter confidence levels of least squares solutions (see, for 
instance [Press et aT]|l992l) . The method consists in generating 
many different data samples by random resampling with repe- 
titions (bootstrapping) the available data, performing the min- 
imization procedure for each sample and deriving confidence 
intervals from the resulting distribution of parameters. The ad- 
vantage, with respect to the errors derived from least squares 
solutions, is that a plot of the parameter distributions directly 
shows the effect of inter-parametric coiTelations, and that the 
confidence intervals are not linked to Gaussian distributions of 
the residuals. In the parameter space the one-sigma confidence 
levels, for instance of a pair of correlated parameters, is defined 
as the projection on the parameter axes of the contour containing 
68.3 % of the bootstrap solutions. 

In this paper we used bootstrap resampling to estimate the 
errors on radial velocities, as obtained from profile fit. This is 
a standard application of the technique: the fitting procedure is 
repeated (1000 times), after bootstrapping the fit residuals. In 
principle the same technique can be used also to estimate the 
uncertainties of the light and radial velocity curve fit, with the 
advantage of deriving more realistic eiTors than that the for- 
mal ones and of gaining insight on parameter coiTelations. It is 
not straightforward, however, to apply the bootstrap technique 
to the complex light curve solution process, because a com- 
plete minimization process should be performed for thousands 
of b ootstrap data-sets. A simpler approach was suggested by 
iMaceroni & Rucinskil (11997 ) who applied bootstrap resampling 
within the minimum already established by a single iterated so- 
lution (that is, using only one set of residuals and parameter 
derivatives). This approach, as discussed in detail in the above- 
mentioned paper, provides underestimated uncertainties with re- 
spect to a full bootstrap procedure, as it is applied already in a 
minimum, but nevertheless, the confidence level we estimated 
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Fig. 7. The contour enclosing 68.3% of all bootstrap solutions 
(1-sigma uncertainties) for two pairs of parameters: upper panel 
/ - Qi, lower panel Qi - Teff ^- The solid cross is the median 
of bootstrap solution values; the arm length is equal to the un- 
certainty on the single parameters, not taking correlation into 
account (i.e the interval of the parameter containing 68.3 % 
of solutions). The broken-line cross is our final solution with 
arms equal to the standard errors. The small offset is due to the 
fact that bootstrap solutions, because of repetitions, are obtained 
from a curve with a diff'erent (and worse) sampling than the orig- 
inal light-curve. The projections of the contours on the axes pro- 
vide the uncertainties of (correlated) parameters mentioned in 
the text. 



for the parameters are always larger than the errors from least 
square minimizations and have a clear statistical meaning. 

The bootstrap uncertainty intervals were, therefore, obtained 
repeating, in this case, 10000 times the last iteration of the fit, 
after random permutation of the residuals. Fig. |7]shows two 2-D 
projections of the n-dimensional distribution of the resulting pa- 
rameter sets (i - Q.\ and f2i - Teff.a)- The figure allows a direct 
comparison between the standard deviations from the PHOEBE 
last iteration and two confidence intervals derived from boot- 
strap: that of each single parameter, not taking coiTelation into 
account, i.e. its interval containing 68.3% of solutions, and the 1- 
sigma confidence level of parameter pairs (the contours), clearly 
showing the effect of correlations (see the plot for the pair Qi 
- Teff,2)- The projection of the contours on the axes provides a 
measure of the uncertainty taking parameter correlations into ac- 
count. 
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The error estimates obtained by this simplified bootstrap pro- 
cedure are somewhat smaller but still the same order than those 
from the heuristic scanning of next section. Taking the largest 
of the full width 1-sigma intervals (as bootstrap errors are not 
symmetric) among various projections and comparing with the 
uncertainties of the heuristic scan of Table |4] we have, 0.02 vs 
0.06 for / (bootstrap vs heuristic scan), 0.01 vs. 0.016 for 
0.01 vs. 0.04 for ^2 and, 40 vs 100 for T^ff^2- The formal errors 
of the fit are typically several times smaller 

The upper panel of Fig. [T] suggests as well a complex shape 
of the cost function along inclination, as the 1 -sigma contour cor- 
responds to a bimodal distribution of the solutions. That feature 
is fuUy revealed by the heuristic scan. 

3.3. Heuristic scanning of parameter space 

Heuristic scanning (iPrsa & Zwitterl l2005h is a Monte-Carlo 
based method to map parameter correlations and cost function 
degeneracy. Scanning produces an A^-dimensional cost function 
map, where is the number of adjusted parameters. Projections 
in that A^-dimensional parameter space reveal solution degener- 
acy. 

We performed the scanning in two sequential steps: the first 
one along inclination, fitting Qi, Q.2, Teff, 2 and computing Li, 
xld and y^o, and the second one along the primary surface po- 
tential, fitting and computing the same parameters except for the 
inclination. 

Because of the unprecedented accuracy of CoRoT light 
curve, a typical fii degeneracy with inclination, surface potential 
Q2 and effective temperature Tsff 2 is essentially broken. There is 
a clear minimum both in inclination (Fig. [8] top left panel) 
and primary surface potential (top right panel) that narrowly 
constrain the uniqueness of the solution. It was very instructive 
to observe the behavior of inclination with the mass ratio: the 
cost function profile varies quite substantially with the changing 
value of q. Although the location of the global minimum remains 
unaffected, adjacent local minima gain on depth, predominantly 
at the expense of "sacrificing" the fit to the secondary minimum 
to better match the abundant out-of-eclipse regions. This indi- 
cates that the model we are using to fit the data is pushed to its 
limit and systematic errors begin to have a significant impact on 
the overall error budget of the solution. 

4. Physical properties of HD 174884 

4.1. Comparison w/fh evoiutionary modeis 

The comparison between the physical parameters obtained for 
HD 174884 components and the stella r models was done ac- 
cording to the method d escribed in lMiglio & MontalbanI (120051) 
and iMiglio et al ] (l2007h . That consists i n applying a Levenberg- 
Marq uardt gradient descent algorithm (iBevington & RobinsonI 



12003 ') to minimize a quality function which describes the dis- 
crepancy between the models and the observables: 



No 



X 



iO' 



,obs 



0l}ieoyl 



1=0 



(2) 



The procedure iteratively adjusts free parameters of the models, 
yielding O'f^™, until best fit is reached. In our case the "observ- 
ables" {O"'"') are the system parameters derived in the previous 
sections: ^j, Teft- 1, the effective temperature difference be- 
tween components, {IsJ^s), the mass ratio (q), the radius ratio {k) 
and metallicity. The corresponding errors are cr"*'. Though we 



initially chose as observables the absolute values (masses, radii, 
temperatures), we finally preferred to use the parameters of the 
primary component and their ratios or differences with respect to 
the secondary, since the light curve solution is mainly sensitive 
to relative values. For the uncertainty on AFeff we assumed, on 
the basis of the light curve fit results, a value of of 100 K. As the 
available data were not accurate enough to derive metallicity, we 
assumed different values of [M/H] with an uncertainty of 99%. 

The theoretical values (O'f"'") are obtained from stel- 
lar evolution modeling with the code CLES (Code Liegeois 
d' Evolution Stellaire, IScuflaire et ani2008bl) . In all model com- 
putations we used the m ixing-length theory (MLT) of convec- 
tion ( Bo hm-VitensellT95 8) and the most rece nt equ ation of state 
from OPAL (OPAL05, Rogers & Navfono vll2002l) . Opacity ta- 
bles are those from OPAL (Iglesi as & Roge rs'1996') for two dif- 
ferent solar mixtures, the standard one from Grevesse & Noeli 
( 19 9ji GN93) and the recently revised solar mixture from 
Aspiund et al.1 (I2005L AGS05) . In the former (Z/X)q=0.0245, 
in the latter (Z/Z )p)-0.0167. The s e tabl es are extended at low 
temperatures with iFerpuson et alj ( l2005h opacity values for the 
correspondi ng metal mixtu res. The surface boundary conditions 
are given bv lKuruc^ (Il998h atmosphere models. 

The parameters of the stellar model are: mass, initial hydro- 
gen (Xi) and metal (Zi) mass fractions, age and two convection 
parameters ((Tmlt and the overshooting parameter OV). All of 
them, or just a subset, can be adjusted in the minimization. The 
value of ffiMLT was kept fixed, adopting the solar value of 1.8, 
as it has no relevant effect on the evolutionary tracks of models 
with masses in the domain of interest here. 

Different minimizations were done with X; as a free pa- 
rameter or fixed to different values between 0.69 and 0.74. 
Similarly, the overshooting parameter, which describes the size 
of the extra-mixed region close to the convective core, was kept 
fixed at values between 0.0 and 0.4, but in some cases left free. 
Binarity puts an additional constraint on chemical composition 
and age, which can be assumed to be the same for the com- 
ponents. Therefore, the number of free parameters varied from 
four i^i, Z, age) to six : ^2, X, Z, OV and age. We 
performed two sets of minimizations, the first with the same 
physical description of components (essentially same value of 
overshooting parameter), and the second allowing a different 
parametrization of extra-mixing. 

Regardless of these physical details all computations show 
the same trend: if the Areft is well fitted (within 1-cr) the mass 
ratio is larger than the observationally derived value and the de- 
viation is typically between 2 and 3cr . This fact can be explained 
considering that for two stars with the same initial chemical 
composition and physical processes, fixing the effective temper- 
ature difference at given age is essentially equivalent to fixing the 
mass difference. Due to correlation, any combination of AT" and 
q yields estimates on q that are inconsistent with the light and 
radial velocity curve solution. This behavior is evident in Fig.|9] 
which shows the iso-^^ contours in the q-AT^ff plane. This plot 
was built using the several thousands binary models computed 
in the minimization processes. If the constraint on AFeff is re- 
laxed, the fit within Icr for all the other constraints is achieved, 
but Areft is ~ 1600 K (Model A of Table IS. Among the mini- 
mizations performed with same physics, AFgft - 1310 K is the 
minimum difference of effective temperature that we were able 
to achieve, keeping the mass and radius ratios within 2 cr (1.6cr 
for q and 1.3 cr for k) (Model Bl). These best fits were obtained 
from computations without overshooting. Values of OV larger 
than zero generally lead to smaller primary masses and hence 
larger q values. 
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Fig. 8. The results of heuristic scanning. Top two panels depict the projection of the cost function space along inclination (left) 
and primary surface potential (right). The bottom two panels depict correlations between the secondary surface potential (left) and 
secondary effective temperature (right). A gray stripe encompasses 1% deviations from the minimum cost function value in the 
primary surface potential, obtained by fitting a parabola to the top right panel. Solid lines represent the best-fit solution, while 
dashed lines represent parameter hmits as determined by the plausible primary surface potential values. 

Table 5. Best fit models. The values of the "observables" , O?**, are given for comparison in the note. 



Model 


^\ 


^2 


Z 


X 


Age 


OV2 


^1 


^2 




9 


k 






(Mo) 


(Mo) 






(Myr) 




(Ro) 


(Ro) 


(K) 






(K) 


A 


4.05±0.09 


2.74±0.06 


0.014±0.003 


0.728+0.090 


121+49 


0.000 


3.76 


2.03 


13189 


0.676 


0.541 


1631 


AC 


4.04+0.11 


2.75±0.08 


0.011+0.003 


0.720* 


112+9 


0.062 


3.68 


1.97 


13737 


0.682 


0.537 


1520 


Bl 


3.99±0.08 


2.74±0.06 


0.011+0.002 


0.740* 


131+9 


0.000 


3.77 


2.00 


13097 


0.688 


0.530 


1310 


B2 


4.02±0.08 


2.77±0.06 


0.010+0.002 


0.720* 


115+8 


0.000 


3.77 


2.00 


13614 


0.689 


0.530 


1338 


C6 


4.01±0.08 


2.76±0.06 


0.012+0.002 


0.720* 


118+9 


0.200 


3.77 


2.00 


13339 


0.688 


0.530 


1295 


Cl 


4.05±0.11 


2.78±0.08 


0.018+0.010 


0.738+10.2 


131+57 


0.654 


3.77 


2.02 


12668 


0.685 


0.536 


1305 


C4 


4.04±0.11 


2.77+0.08 


0.015+0.007 


0.740* 


130+13 


0.355 


3.77 


2.02 


12870 


0.686 


0.535 


1320 


C2 


4.00±0.08 


2.75±0.06 


0.013±0.002 


0.740* 


131+10 


0.200 


3.77 


2.00 


12978 


0.687 


0.530 


1277 


C5 


4.00±0.08 


2.75±0.06 


0.012+0.002 


0.740* 


131+1 


0.100 


3.77 


2.00 


13045 


0.688 


0.530 


1293 


C7 


3.99+0.08 


2.75±0.06 


0.011+0.002 


0.740* 


130+95 


0.050 


3.77 


2.00 


13078 


0.688 


0.530 


1298 


Of ': J- 


=4.04 + 0. 11 =3.77 


+ 0.03;reffi = 


13140 + 1500; q 


= 0.674 


+ 0.009; k 


= 0.541 


+ 0.013; Areft 


= 1096 + 100 (or 1240+10 



see text). 

For all models OVi = 0. Model AC belongs both to group A and C. Several C-type models are listed to illustrate the effect of adjusting OV2. 
*' Fixed value. 



Is it possible to reconcile the observed values of q and ATeff 
with those from the theoretical models? To answer this ques- 
tion we considered the possible alternatives and the weak points 
of comparison with the models. An evident inconsistency be- 
tween the stars and the corresponding models is that CLES com- 
putes spherical non-rotating models, while the primary and sec- 



ondary component of HD 174884 rotate, respectively, at 25% and 
12% of their break-up critical velocity, and are significantly non- 
spherical. The distortion of the primary {Rpoie ~ 2.5% smaller 
than 7?eq) is mainly due to rotation. It is well known that at given 
mass, rotating stars appear to be cooler than non rotating ones, 
or - expressed in different terms - at given elfective temperature. 
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non-rotating models will assign lower stellar masses than rotat- 
ing ones. The difference of effective temperature between rotat- 
ing and non-rotating models depends on the rotational velocity, 
inclination of t he rotation axis and che mical composition of the 
star (see, e.g., iMaeder & Mevnete004i) . An estimate of the av- 
erage effect of rotation on the ob served effective temperatur e for 
moderate rotators is provided bv'Bas tian & de MinkI ( |2009|) and 
their results are in agreement with those of the Geneva group 
for a 4 M0 star (Meynet, private communication). According to 
the above-mentioned authors, the change in T^ff of a rotating star 
with respect to a non-rotating one can be expressed as : 



T,«(ji)/T,b(0) = 1 - ar]\ (3) 

where 77 is the fraction of the critical break-up rotation rate, 
TeffiO) is the effective temperature of the non-rotating model, 
and fl is a coefficient which changes from 0.17 at zero age to 
0.19 at the end of the main sequence. In the case of HD174884, 
the apparent decrease in temperature due to rotation amounts to 
~ 160 K for the primary and only 20 K for the secondary compo- 
nent. As consequence, the observed effective temperature differ- 
ence should be smaller than that of non-rotating models, and that 
could explain the higher AT^ff values from best fit models with 
respect to that from the light curve solution. In this hypothesis, 
the constraint on AT^ff to be used with our non-rotating models 
shall be increased from 1096 K to ~ 1240 K. This value, how- 
ever, is still somewhat lower than that obtained in our best fit. 

The value of the overshooting parameter, its dependence on 
the stellar mass, and the physical origin of the extra-mixing are 
still a matter of debate. A different value for the two components 
cannot be completely excluded and the hypothesis OV] - OV2 
could be dropped. We tried, therefore, to satisfy our observa- 
tional constraints by using fixed and different values of OVi and 
0V2- We also performed computations fixing OV for the pri- 
mary and varying the OV for the secondary component, and - 
as well - keeping fixed OV2 and deriving OVi. The correspond- 
ing models are labeled as models C in Table |5] An overshooting 
parameter higher for the secondary than for the primary slightly 
improves the fit with respect to model B. From a physical point 
of view, we could justify a secondary with a stronger mixing 
in the center, even if its rotational velocity is lower than that of 
the companion, by invoking an important gradient of the inter- 
nal rotational profile produced by braking in the synchronization 
process. However, the improvement in the fit is not significant 
enough to justify the introduction of an additional parameter. 

The age of the system is similar for all models and has an 
average value of 125 + 7 Myr 

Finally, we can ask ourselves if our conclusions could de- 
pend on the evolutionary code we used for computing stellar 
models. A convincing answer to this question can be found in 
the extensive study which was caiTied out by the CoRoT/ESTA 
(Evolution and Seismic Tools Activity) group, in view of the 
CoRoT mission and the foreseen analysis of seismic data of 
the highest ever achieved accuracy. The team made an accurate 
comparison among different evoluti onary and pulsation c odes 
for selected test cases. Their results (iLebreton et al.ll2008l and 
other contributions appearing in the same special issue) allow 
us to conclude that in the part of the parameter space of interest 
here (MS stars of intermediate mass) stellar evolutionary codes 
adopting the same physics yield the same fundamental parame- 
ters and the same internal structure. 
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4.2. Orbital evolution 

The position of HD 174884 in the period-eccentricity diagram of 
close binaries is shown in Fig. [10] The plot is based on the cat- 
alog of Hegedus et al. (2005), with the addition of HD 313926, 
a short period (P=2?27) high eccentricity (e=0.21) eclipsing bi- 
nary formed by two B-type stars, which was rec ently discov- 
ered by the MOST satellite dRucinski et alj|2007l) . The major- 
ity of B-binaries in the figure have main sequence components. 
HD 174884 has the second highest eccentricity among binaries 
wit h B type components. 

iRucinski et al.l (l2007h suggest that the upper envelope for 
B type stars in the period-eccentricity plot might be flatter than 
that of the later stellar types. This can be understood in terms 
of stronger dissipation, and hence more efficient circulariza- 
tion mechanisms, for s tars w i th larg er fractional radii, see for 
instance the review of (l200l . That is at the origin of 

the overall trend in the plot (lower eccentricity for shorter pe- 
riods), besides, at fixed period, less massive (unevolved) stars 
have smaller fractional radius and hence longer circulariza- 
tion timescales. The evidence of the second effect is, however, 
marginal in our updated plot, as HD 174884 occupies a region 
previously empty of B-stars. 

The relation between eccentricity and the fract ional radius of 
early type binaries has been extensively studied. iGiuricin et alJ 
( 1984) analysed a sample of galactic binary and North & Zah^ 
(2003) more homogeneous ones in the Magellanic Clouds. All 
samples indicate the existence of a "cut-off fractional radius, 
above which all systems have circular orbits. This value - prac- 
tically independent of chemical composition - is ^ 0.25. The 
fractional radius of HD 174884's primary is ri ^ 0.2, i.e. a value 
consistent with its ellipti c orbit. Its posit i on in t he eccentricity vs. 
fractional radius plot of IGiuricin et al.l (1 19841) (not shown here 
for the sake of brevity) marks again the upper envelope of the 
galactic sample for that period. 

The same mechanism producing orbit circularization is at 
the origin of spin-orbit synchronization on timescales typically 
a few orders of magnitude shorter than that for circularization, 
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Fig. 9. The iso-;^^^ contours in the q-Arefr plane for the 
Levenberg-Marquard minimizations of the deviations with re- 
spect to theoretical models. Lower values correspond to darker 
colors. The plot clearly shows the correlation between the two 
parameters. The is computed with respect to the AFobs value 
obtained after coiTection for the effect of rotation (see text). The 
original value is also shown, as a staiTed dot. The broken line 
boxes enclose 1- and 2-cr intervals. 
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P(d) 



Fig. 10. The period - eccer itricity diagrarn , for cl ose binaries ac- 
cording to the catalog of iHegedus et al.l ("2005). Filled circles 
correspond to B stars, open circles to A-G stars. The star sym- 
bol indicates the location of HD 174884, the triangle that of 
HD 313926 dRucinski et al.ll2007l see text). 
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Fig. 11. Sections of the original light curve, phased according to 
the orbital period and vertically shifted for better clarity. The 
seven observed orbital periods go bottom-up. The oscillatory 
pattern of the residuals is clearly seen in each curve. 



as - in presence of efficient angular momentum transfer - the 
ratio of the two timescales equals that of rotational and orbital 
angular momentum content. As long as the orbit is elliptical, 
the component rotation will synchronize with the orbital motion 
at periastron. Since v sin / was estimated from spectral analysis 
(visin/ - 109km s"',V2 sin/ = 60km s"') and the eccentric- 
ity is well constrained from spectrophotometric modeling, we 
can readily verify whether this is the case. A theoretical value 
dHutll 19811) of the synchronicity parameter for synchronous ro- 
tation at periastron is F - Vo'^'^'ejTo'^-ej^, and for our system 
F - 1.92, while spectroscopic analysis yields F[ = 2.0 + 0.2 and 
F2 - 2.3 + 0.2. Both stars thus seem to be rotating marginally 
super-synchronously. We can check if the dynamical state of 
the system is in agreement with the expected timescales. If 
we follow Zahn's formalis m (see the above-mentioned review 
and IClaret & Cunhal 1 19971) to estimate the circularization and 
synchronization timescales of HD 174884, a dopting a value of 
B = 0-2 for the fractional gyration radius (IClaret & Gimenea 
l2005b and log £2 — -6.6 for the value of the second torque con- 
stant ( IClaret & Cunhal 1 19971) . we get: Tdrc = 8.9 • 10^ yr ( or 
35% shorter considering the contribution of both components), 
Tsynci - '!"drc/120 = 7.4 ■ 10^ yr and Tsync,2 - 4Tsync,i. While the 
relations providing these estimates are derived under the hypoth- 
esis of small deviation from circular orbit and synchronism, the 
resulting figures are indeed in good agreement with the dynami- 
cal status of our system. 

The fact that in the same period interval of the e-P plot both 
circul arized and elliptic orbits can be found is interpreted (Zahn 
l2005h as due to m uch higher efficiency of tidal damping when 
resonance locking (IWitte & Savoniielll999i) takes place. As this 
event is very sensitive to the binary parameters it might appear 
only in some systems among those of similar period, depending 
on the other stellar parameters. According to this interpretation 
the orbital evolution of HD 174884 has not been driven by reso- 
nance locking. 



4.3. Pulsational properties: analysis of iigint-curve fit 
residuais 

The residuals from the phased light curve fit, appearing in the 
lower panel of Fig.|2] clearly show a phase-locked pattern, which 
can be described in terms of two components: a complex irregu- 
lar structure at phases close to (and within) primary eclipse and 
a multi-periodic oscillatory behaviour at other phases. Their dif- 
ferent nature is suggested by the fact that only the latter is present 
in the original light curve (see Fig.lTTIl. 

We interpret the first component as due to systematic effects 
caused by the employed EB model, e.g. to small deviations in 
model atmospheres (and consequently in the limb darkening de- 
scription), and to the description of stellar surfaces by a finite, 
though large, number of elements. These effects are expected 
to be visible during eclipses, when one stellar disc eclipses the 
other Those related to the surface discretization are typically be- 
low 0. 1 mmag, but can be larger at some critical phases. 

We attempted to identify those segments of the code that 
need to be restated in order to increase the model accuracy. 
An important one is the above-mentioned surface discretization, 
which currently relies on equidistant partitioning of the surface 
along stellar latitude and longitude. This introduces a "seam" 
along (fi = Q where all surface elements are perfectly aligned. We 
have been devising a modified approach where surface elements 
are determined by equidistant partitioning along the length of 
the equipotential and shifted by half the element width at = 0. 
Preliminary tests indicate that the systematic discrepancy in the 
vicinity of the primary minimum residuals (around tp - -0.025) 
is greatly reduced when employing the new strategy. However, 
significant testing is required to identify any potential problems 
with the new scheme. 

The second component of the residuals, the oscillatory pat- 
tern which can be spotted by eye in the original curve, is cer- 
tainly model-independent and is due either to the intrinsic vari- 
ability or to an instrumental artifact. The latter hypothesis is, 
however, unlikely, given the perfect phasing with the binary 
period over more than seven cycles. Furthermore, none of the 
known frequencies of instrumental origin is close to those found 
in the pattern components. To extract the individual frequencies 
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Table 6. Results of Fourier analysis 



Fn 


Frequency 


Amplitude 


Phase 


S/N 


Remark 




(d ') 


(10-3) 








Fl 


0.273444 


0.452 ±0.002" 


0.3848 ±0.0006 




/o,b, fixed 


F2 


0.5470 +0.0002 


0.284 ±0.002~ 


0.840 ±0.001 


14.2 


- 2/orb 


F3 


2.1858 +0.0004 


0.120 ±0.002 


0.349 ±0.002 


12.0 


- 8/orb 


F4 


3.5534 +0.0004 


0.111 ±0.002 


0.822 ±0.003 


11.2 


- 1 3/orb 


F5 


1.0894 ±0.0004 


0.097 ±0.002 


0.840 ±0.003 


6.2 


- 4/orb 


F6 


0.0631 ±0.0004 


0.116 ±0.002 


0.931 ±0.003 


7.2 


long term trend 


F7 


0.8206 ±0.0005 


0.091 ±0.002 


0.570 ±0.003 


7.8 


- 3/orb 



" all uncertainties are the fit formal errors from Period04 



we pe rformed a Fourier analysis with Period04 (iLenz & Bregerl 
I2005h . For this purpose a complete set of residuals, in the time 
domain, was obtained by subtracting the fitted model from the 
~ 60000 points of the original light curve. The analysis is made 
more complex by the presence in the time series of residuals pro- 
duced by the above-mentioned systematic effects, which gener- 
ate power at the orbital frequency (/orb = 0.273444 d"') and its 
overtones. The analysis was therefore performed on two differ- 
ent time strings: the complete set of residuals and a "cleaned" 
sub-sample, obtained by deleting the points the phase interval 
-0.1 < (f < 0.05. Given the high total number of points it was 
possible as well to analyze only some sections of the light curve 
residuals, e.g. each individual period. 

All the analyses provide essentially the same result: a clear 
detection of two frequencies, corresponding to 8 and 13 times 
the orbital frequency and marginal detections of other orbital 
overtones. Table|6]and Fig [12] summarize the results of the anal- 
ysis for the "cleaned" set of residuals. The first frequency cor- 
responds to /orb and was kept fixed. The first two frequencies 
are actually needed to fit the slow trend, visible in the lower 
panel of Fig.|2l due to the systematics introduced by subtraction 
of the light curve model. Similarly the low frequency F6 fits a 
long term residual trend, a left-over of the detrending procedure 
which was applied to the original curve. Frequencies F3 and F4 
correspond to 8 /orb and 1 3/orb- 




S 6.0x10"= 
< 
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2.0x10 



f (c/d) 

Fig. 12. The Fourier spectrum of the "cleaned" residuals after 
pre-withening of the first two frequencies in Table |6l The dotted 
lines mark the two frequencies 8 /orb and 13 /orb- The dashed line 
is the spectrum of noise, computed by Period04 on the residuals. 
The inset shows the spectral window. 



The detections are very clear, reaching amplitudes over ten 
times the noise level. More multiples of the orbital frequency 
are present, but at lower amplitude so we consider them to be 
only marginal detections. Intrinsic variability at multiples of the 
orbital frequency are a phenomenon pointing towards tidally in- 
duced oscillations or free oscillations affected by the close bina- 
rity. 

Even though several pulsating stars in close binaries have 
been found (e.g. lAerts & Harmane3l2004l) . oscillation frequen- 
cies at multiples of the orbital frequency have only seldomly 
been detected. The most likely explanation is that they are due 
to resonant dynamic tides which can lead to detectable ampli- 
tudes. As discussed by, e.g., Wil lems & AertsI (j2002), the shape 
of tidally induced observables (light curves or radial velocities) 
can vary from very iiTegular for orbital periods away from a res- 
onances with free oscillation modes to sinusoidal for orbital peri- 
ods close to a resonance with a free oscillation mode. It is some- 
what surprising that we find clear multiples of orders 1 3 and 8 for 
HD 174884, while the lower order multiples are not so markedly 
found. Interestingly, such a case also occurs for the slowly pul- 
sating B star HD 177863, where an oscillation frequency of ex- 
actly ten times the orbital fre quency was found and interpreted 
as due to resonant excitation dDe Cat et al.ll200(il) . 

According to theory, tides can efficiently excite the free os- 
cillation modes of the star which are close to the tidal frequency 
(and its multiples in eccentric orbits see, e.g., Wil lemsl (20031) . 
The dominant tidal term is associated with the spherical harmon- 
ics of degree f = 2. Besides, given the typical values of the tidal 
frequencies in close systems, the stellar free oscillation modes 
should be g -modes, and typically of high radial order, n > 10 
(IZahnll2005h . 

By means of th e code LOSC (Liege OScillation Code 
IScuflaire et al.ll2008ah we computed the eigenfrequencies of the 
primary star for all the models within 2-cr in mass and radius and 
^Tgff < 1600 K. The purpose was to check the properties of the 
{ = 2 modes close to the F values of Tabled LOSC includes the 
effects of rotation on the frequencies in a perturbative approxi- 
mation, and the correction terms are limited to the first order in 
Wrot (e-g- Ledoux 1951)- This approach is acceptable when the 
effects of the centrifugal force and of Coriolis acceleration are 
small. The former is estimated by the ratio (e) between the rota- 
tional period and the dynamical time , the latter by the ratio (ju) of 
rotation to oscillation frequencies. For the HD 174884 primary 
component e ~ 4 x 10'^ and = 0.25 for F3, and fi = 0.15 for 
F4. As shown by iMigUo et"al] (I2008h . the effect of rotation on 
the internal chemical composition profiles can also lead to slight 
changes of the eigenfrequencies. 

Since our aim is not to fit the observed frequencies, but to 
verify that their value is well in the range of f - 2 g-modes 
with radial order n of ~10 or larger, we consider the first order 
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approach sufficient. The comparison between observed and the- 
oretical frequencies indicates, indeed, that the lower frequency 
F3 corresponds, for the models described in Table |5] to n values 
between 8 and 10, and the higher one to slightly larger values. 
We can conclude, therefore, that the frequencies we find in the 
residuals are compatible with tidally forced oscillations. 

As a further check of the light curve solution robustness, we 
pre-withened the original curve by the two frequencies F3 and 
F4 and recomputed the best fit. The new solution is essentially 
the same as the old one: the largest difference between best-fit 
parameters before and after pre-withening is 0.4%. The sub- 
traction of the harmonic pattern does not substantially improve 
the correlation of residuals, because their dominant component, 
which we interpreted as due to model systematics, is essentially 
left unchanged by the procedure. 

5. Summary and conclusions 

This work has been carried out with two aims in mind: to achieve 
a thorough description of an interesting binary and to test the 
performance of current binary star modeling on data of unprece- 
dented accuracy, as those available for the CoRoT seismo-field 
targets. HD 174884 lends itself to both purposes: its peculiar 
light curve suggests an unusual system configuration, worth a 
detailed study, and it is as well quite stable over the seven ob- 
served orbital periods so that, if possible systematic errors aris- 
ing from binary modeling were present, they would not be hid- 
den by transient intrinsic phenomena. (The latter is indeed the 
case of the other fu lly analyzed binary seismo-target of CoRoT 
first runs, AU Mon (iDesmet et al.ll2009b whose study reveals the 
presence of variable circumbinary material). 

Given the presence of a grazing eclipse, we were concerned 
about solution uniqueness, a problem which is sometimes over- 
looked, and we put a great effort in getting a sound estimate of 
parameter uncertainties, essential for sensible comparison with 
evolutionary stellar models. 

We can summarize our main results as follows: 

- due to the high accuracy of CoRoT data, we were able to un- 
ambiguously derive parameters, such as inclination and com- 
ponent sizes, which would otherwise be poorly constrained 
in a system with grazing eclipses. This fact, the excellent 
agreement between the results of the different methods ap- 
plied for the analysis of the spectra, and the extensive check 
of solution uniqueness allowed a sound estimation of the 
physical parameters of HD 174884 and their uncertainties. 

- The comparison with stellar models is quite satisfactory 
but for the temperature difference between the components, 
which, for the best fitting models, is typically a few hundred 
degrees higher than that from the light curve solution. We 
interpreted that as due to the comparison with non rotating 
stellar models. We could as well improve the agreement with 
models introducing differences in the physical processes act- 
ing in the components, such as overshooting, but Occam's 
razor arguments favor the simpler, though not exhaustive, ex- 
planation. Increasing the free parameters yields, in fact, only 
a marginal improvement. The comparison with theoretical 
models provided as well an estimation of the system age and 
some indication on the component chemical composition. 

- The dynamical properties of HD 174884 are in very good 
agreement with the predictions of Zahn's theory of circular- 
ization and spin-orbit synchronization. The high value of the 
eccentricity for the period suggests that resonance locking 
has not been at work in the system. 



- A few frequencies are clearly detected in the residuals of the 
light curve fit at multiples of the orbital frequency. We ten- 
tatively interpret them as resonantly excited pulsations. This 
hypothesis is strengthened by the fact that many g-modes of 
high radial order and degree ( - 2 exist in the frequency 
range of tidal frequencies. Among the very few detections of 
tidally induced pulsations in binaries, HD 174884 is the only 
case with two well characterized early type MS components. 

Finally, it is evident from our results that when dealing with 
observation accuracy of a few tenths of mmag we are close to the 
level at which systematics from the model has a relevant impact 
on the light curve solution. This indicates that the models shall 
be updated, especially in view of the still higher accuracy curves 
which will be obtained by curre nt and future space missions such 
as Kepler dBorucki et al.l 120081) . Similarly, the physical proper- 
ties derived by data of such accuracy should be compared with 
theoretical models including stellar rotation. 
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